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The presence of substructures in dark matter haloes is an unavoidable consequence of the cold dark 
matter paradigm. Indirect signals from these objects have been extensively searched for with cosmic 
rays and 7 rays. At first sight, Cherenkov telescopes seem not very well suited for such searches, 
due to their small fields of view and the random nature of the possible dark matter substructure 
positions in the sky. However, with long enough exposure and an adequate observation strategy, 
the very good sensitivity of this experimental technique allows us to constrain particle dark matter 
models. We confront here the sensitivity map of the HESS experiment built out of their Galactic scan 
survey to the state-of-the-art cosmological N-body simulation Via Lactea II. We obtain competitive 
constraints on the annihilation cross section, at the level of 10 -24 — 10~ 23 cm 3 s _1 . The results are 
extrapolated to the future Cherenkov Telescope Array, in the cases of a Galactic plane survey and 
of an even wider extragalactic survey. In the latter case, it is shown that the sensitivity of the 
Cherenkov Telescope Array will be sufficient to reach the most natural particle dark matter models. 

PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,95.30.Cq 



I. INTRODUCTION 

In the current cosmological paradigm, cold dark mat- 
ter (CDM) dominates structure formation. The haloes of 
galaxies and clusters of galaxies are assembled through 
the merging of a huge number of smaller structures. Most 
mergers are incomplete and large CDM haloes, e.g. the 
one around the Milky Way, harbor an enormous popula- 
tion of subhaloes, which are a record of its assembly his- 
tory. Some particle physics models beyond the standard 
model predict the existence of these new weakly interac- 
ting massive particles (WIMPs) that could well form the 
cosmological dark matter (DM). Should these particles be 
thermally produced in the early Universe, their comoving 
density would have been regulated by self-annihilations. 
The freeze-out of this reaction in the early Universe leads 
to the establishment of the nowadays observed CDM den- 
sity. In dense enough structures, the very same annihi- 
lation process can efficiently convert particle DM mass 
energy into high energy standard model particles. The 
CDM structures represent at least one-tenth of the total 
halo mass and are privileged targets of searches for DM 
particle annihilations. The dense regions worth conside- 
ring are the Galactic center and galactic subhaloes, in- 
cluding a small subset of them which harbors the dwarf 
satellite galaxies of the Milky Way. The principle of indi- 
rect searches for DM through 7 rays is to search for the 
7-ray emission following the hadronization and/or decays 
of these cxotically produced standard model particles. 

Strategies for searching this dim radiation include tar- 
geted searches or wide-field surveys. Thanks to their very 
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large collection area, imaging atmospheric Cherenkov te- 
lescopes (IACT) are very well suited for deep observa- 
tions of selected sources. Constraints on particle DM mo- 
dels have been obtained by the HESS experiment from 
the observation of the Galactic center pj, the Sagit- 
tarius dwarf galaxy Q, and the Canis Major overden- 
sity Q and by MAGIC from the observation of the Per- 
seus galaxy cluster j3] and Milky Way satellites Draco [f| 
and Willman 1 @. The VERITAS experiment obtained 
constraints from the observation of Draco, Ursa Minor, 
Bootes 1, and Willman 1 Q and the Whipple Colla- 
boration got constraints from Draco, Ursa Minor and 
M15 Q. In this paper we focus instead on nontargeted 
searches with Cherenkov telescopes. Space-based instru- 
ments such as the Fermi satellite can much more easily 
perform blind searches for DM subhaloes with a regular 
scanning of the entire sky thanks to their large field of 
view. Some prospect regar ding searches with Fermi have 
been conducted in Refs. |9l— Tlll| and recently the first ana- 
lyses based on actual data were released (Refs. 
In [l4j], HESS data from the Galactic plane survey have 
been used to perform for the first time a blind search for 
DM substructures with a wide-field survey with IACTs 
(the structures were DM spikes around intermediate- 
mass black holes ). However, this substructure scenario is 
rather optimistic since the abundance and the properties 
of intermediate-mass black holes and of the DM spikes 
around them remain practically unconstrained. Here we 
investigate constraints derived using the HESS Galactic 
survey and the conventional CDM subhalo distribution 
obtained by the cosmological N-body simulation Via Lac- 
tea II (VL-II) [la ]. As we shall see in the following, the 
outcome of this study is based only on the numerically 
well resolved distribution of CDM structures in the Ga- 
lactic halo, it is thus quite robust. Also, it does not rely 
on further density enhancements such as e.g. the pos- 
sible formation of intcrmcdiatc-mass black holes. This 
results in weaker but safer constraints on the particle 
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physics models than those of [l4|. To go beyond these 
constraints, it is interesting to consider extended surveys 
that will certainly be performed by the next generation of 
IACT such as the Cherenkov Telescope Array (CTA) [H 
as high priority observations. We thus extend the study 
by extrapolating the current constraints to the sensitivity 
of CTA, the future observatory in this range of energy. 

The paper is structured as follows. In Sec.|II]we present 
the results of the VL-II simulation regarding the subhalo 
population which we are interested in. Section IIIII is de- 
voted to the description of the HESS sensitivity map and 
how the extrapolation to the CTA is done. The results 
inferred from the HESS survey are presented in Sec. IIVI 
Finally, Sec. [V] presents the prospects with the CTA and 
Sec. IVIlis devoted to the conclusion. 



PREDICTIONS FOR THE SUBHALO 
CONTENT OF THE MILKY WAY 



II. 



The CDM subhalo distribution is taken directly from 
the VL-II simulation [15j . one of the largest, most ac- 
curate cosmological N-body simulations of the Galactic 
CDM halo. The particle mass of 4100 Mq allows us to 
resolve small subhaloes (> 10 5 M Q ) throughout the Ga- 
lactic halo and even as close to the Galactic center as 
the solar neighborhood. Much smaller CDM subhaloes 
are expected to survive as well [13], but we will show 
now that for the purposes of this analysis it is not neces- 
sary to include CDM clumps below the VL-II resolution 
limit. The subhaloes resolved in VL-II increase the lumi- 
nosity of the entire haloes as seen by a distant observer 
by a factor of 1.9. Including the smaller unresolved ha- 
loes expected to exist in CDM does increase this factor 
to values between 4 and 15, depending on the abundance 
and properties assumed for these smaller clumps. The 
total luminosity from subhaloes is also increased by si- 
milar factors due to their own subsubstructrue [|, but 
we assume smooth subhaloes throughout this work for 
simplicity. Locally (within 1 kpc of the Solar System) 
this substructure models leads to an enhancement of only 
1.4±0.2 [3. 

The dark matter annihilation luminosity 1 is defined by 
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where p 2 is the squared mass density of the subhalo. The 
smallest clumps, which are still well enough resolved in 
VL-II, have peak circular velocities of V m!iK = 3 kms -1 
and L — 1.7 x 10 5 MqPC -3 , and their local separation 
is (d) = 5.8 kpc; i.e. in a random realization they are 
found at a median distance of Do — 0.5(d) ~ 2.9 kpc 



1. Here we assume a smooth subhalo. Clumps within a subhalo 
do enhance the luminosity but are ignored here. 



from the observer [18j . The 7-ray flux from a given clump 
is calculated from 
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Nj being the integrated number of 7-rays produced in 
a WIMP collision above a given energy threshold E t h, 
m the WIMP mass, and av the velocity- weighted anni- 
hilation cross section. With D = 2.9 kpc, one gets the 
reference flux of $0 = 7 x 10~ 14 cm~ 2 s _1 from such a 
small, but still well resolved, VL-II clump . This value is 
obtained for av = 3 x 10 -26 cm 3 s _1 , m = 500 GcV in 
the case of an annihilation into r pairs. As we shall see 
in the following, the observable flux for HESS is roughly 
$hess = 10" 12 cm~ 2 s~ 1 = 14 x $ and in the best 
case for the CTA is ~ 3 x 10~ 13 cm~ 2 s _1 , which is still 
larger than $o- This means that the smallest clumps that 
the VL-II simulation is able to resolve are already too 
faint to be observed by HESS or the CTA for most ran- 
dom realizations. They would need to fall unusually close 
to the observer to be detected (D < 1.8 kpc), which hap- 
pens in only 18% of all realizations. A HESS detection 
of a smaller CDM clump, below the numerical resolution 
of our simulation is even less likely : In CDM the dis- 
tance D n to the nearest subhalo of above given luminosity 
scales roughly like D n oc L 1 / 3 (see for instance (HI). The 
flux $ oc L/D 2 therefore goes like $ oc L 1 / 3 , i.e. the flux 
from the nearest subhalo in a given bin in logL (or log 
subhalo mass) increases with subhalo mass and we can 
safely ignore small subhaloes below the VL-II numerical 
resolution for the present analysis. 



III. FLUX SENSITIVITY MAPS 

Data from the Galactic plane survey conducted by the 
HESS experiment allowed us to obtain flux sensitivity 
maps to DM annihilations [TH- In R- CI - H3, the 7-ray 
flux sensitivity is calculated in each position depending 
on the acceptance of the detector, the exposure to each 
point of this region of the sky and the predicted 7-ray 
spectrum, according to 
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Here, A g is the effective area for 7 rays, which is a func- 
tion of the gamma energy E 7 , the sky position (b,l), and 
T D bs the observation time at a given sky position. No- 
tice that the flux sensitivity depends explicitly on the 
mass of the DM particle m and on the details of the 
annihilation process. The HESS Galactic plane survey 
allowed us to detect a large population of unidentified 
sources [III H| . None of them exhibit an energy cut-off 
in the covered energy range. Each 7-ray source is best- 
characterized by a power-law spectrum, as expected for 
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Figure 1: Flux sensitivity maps at 90% C.L. for HESS and a CTA-like array. The maps are calculated here for a 500 GeV DM 
particle annihilating with 100% BR into bb. Top : The HESS sensitivity map is calculated with real exposure. Bottom : The 
simulated flux sensitivity map for a CTA-like array is obtained for a flat exposure of lOh. 



7-ray emission from standard acceleration astrophysical 
processes. Consequently it excludes them from being rele- 
vant candidates for DM clumps. No further information 
about the energy spectrum is used because the nonde- 
tection of DM clumps -which is the baseline assumption 
throughout this paper- depends first on the total flux. 
As far as the integrated flux is used, the energy spec- 
trum of the searched sources only influences the results 
indirectly though the choice of the annihilation channel. 
Fig. Q] shows the HESS flux sensitivity map to DM an- 
nihilation for a 500 GeV DM particle annihilating into 
bb. In the following, when other annihilation channels or 
masses are considered, the flux limit is properly resca- 
led. The sensitivity map contains hints of the presence 
of all the discovered sources by HESS. As expected, the 
flux sensitivity decreases at the position of the detected 
sources. The flux sensitivity map can be easily unders- 
tood as follows : if a 7-ray source at a given Galactic 
position gave a larger flux than the value quoted in the 
corresponding bin, it would have been detected. Conse- 
quently, any model predicting a statistically significant 
number of sources of any kind with fluxes above these 
values is excluded. 

An important point concerns diffuse emission, from 
either cosmic rays, the smooth DM halo or unresolved 
clumps. Note that the Cherenkov technique is suitable for 
detecting pointlike sources and slightly extended sources 
or diffuse emission. The reason is that one has to per- 
form an on-off type of background subtraction (see for 
instance [HI). Because the trigger rate can vary from 
one observation to another (depending on the observing 



conditions, the zenith angle, night sky background level, 
etc.), the angular scale one has to keep in mind is the size 
of the effective field of view of an individual run, which is 
about 4° x 4° 2 . Unless a more subtle method is used, the 
off region has to be taken in the same field of view as the 
on region to avoid large systematic errors. This means 
that if a signal is spread over a region that is larger than 
5°, it cannot be observed by Cherenkov telescopes using 
the on-off technique. Although they concern large fields 
of view, the maps of Fig. Q] are built out of a concatena- 
tion of smaller 4° x 4° observations. Concerning diffuse 
emissions on scales below ~ 4°, some has been detected 
by HESS (e.g. in the Galactic center region (22J) and is 
an additional source of background, the point-source sen- 
sitivity within those regions is decreased. This is clearly 
visible when scrutinizing the map of Fig. Q] Any diffuse 
signal that would be pretty much constant all over the 
map -should it reach or even exceed the point-source sen- 
sitivity of the experiment- would actually be erased by 
the background subtraction procedure. In practice it is 
always the case ; for instance the signal from the smooth 
halo and unresolved clumps vary by a factor of at most 
~ 20% on a 4° scale :23| . The only region where it can be 
untrue is the central 1° where the central cusp of the DM 
halo could appear as a slightly extended source, hence 



2. On the CTA, the individual field of view will be larger (> 7°) 
for energies above 100 GeV, but this does not qualitatively affect 
the discussion on diffuse emission 
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producing a large gradient diffuse emission. This is not 
a concern for our study, because of the presence of the 
central astrophysical source. The issue of an annihilation 
signal at the Galactic center has been addressed anyway 
in In a sense, because of the way the sensitivity map 
is built, the issue of large scale diffuse emission is handled 
automatically. This is true as long as the sources that are 
searched for have an extension lower than ~ 2° . As we 
shall see in the following, it is always so for the considered 
case. 

The next generation of IACT will consist of a large 
telescope array (the CTA) (24[. The current effort on its 
design should allow us to improve significantly the glo- 
bal performance of the present generation. The goal is 
to extend the accessible energy range towards both the 
low and the high energies and gain at least a factor of 
10 in flux sensitivity. To extrapolate the current results 
to future observatories, different assumptions are made 
regarding the foreseen CTA characteristics. A conserva- 
tive value of 50 GeV for the energy threshold is assumed. 
The effective area of the instrument is improved by a fac- 
tor of 10 and the hadron rejection by a factor of 2. The 
overall improved capabilities of the CTA will presuma- 
bly allow us to detect a number of new 7-ray sources. 
The construction of the CTA-extrapolated map follows 
the procedure described in [25[. Assuming a supernovae 
remnant source model for the 7-ray emission [26j and a 
radial source distribution [27} , the distance and the 7-ray 
flux are calculated. The Galactic plane is then randomly 
populated according to the spatial distribution of sources 
observed by the HESS survey [H[IJ. So the HESS 7-ray 
flux source distribution is extrapolated to CTA perfor- 
mance. This extrapolation results in the prediction for 
the discovery of a few hundreds of new sources in the 
Galactic survey field of view. The presence of these new 
sources deteriorates the DM flux sensitivity accordingly. 
The resulting projected CTA sensitivity map is shown in 
the lower panel of Fig. [1] for a DM particle mass of 500 
GeV annihilating into 66. A flat exposure of 10 h in each 
position of the map is assumed. This value allows us to 
match the total amount of time for the CTA survey to the 
~ 400 h which were needed by HESS to survey this re- 
gion of the sky. The flux sensitivity for CTA ranges from 
~10 -12 cm" 2 s _1 in the region where the new sources are 
present to a few 10 _13 cm _2 s _1 on average at higher la- 
titudes. 



IV. CURRENT EXCLUSION LIMITS 

The question to be addressed in this section is -for a 
fixed set of particle physics parameters- what the pro- 
bability is for a clump to lie in the survey region with 
a flux larger than the HESS sensitivity at its position 
in the sky. To answer this question, we generate 10 3 dif- 
ferent Monte Carlo realizations of the Milky Way halo. 
Independent realizations are obtained by randomly pla- 
cing the observer at a distance of 8.5 kpc from the Ga- 
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Figure 2: Distribution of the number of subhaloes between 
±3° in Galactic latitude and from -30° to 60° in Galactic 
longitude. The distribution is obtained from 1000 stochastic 
realizations of the Milky Way-like halo from VL-II. The mean 
of the distribution is 168 with a rms of 44. 



lactic center. Over our realizations, the total number of 
subhaloes inside the Galactic survey field of view shown 
in Fig. Q] is 168 ±44, out of the ~ 10 4 resolved subhaloes 
contained in the Milky Way. The distribution obtained 
after these 10 3 virtual experiments is displayed in Fig. [5] 
The probability to find no subhalo in the field of view 
is less than 10 -3 . The spatial distribution of clumps is 
slightly triaxial, with an unknown orientation which is 
not expected to be correlated with the baryonic distribu- 
tion. For that reason, the distribution is wider than what 
is expected from a purely spherical distribution. To scan 
the (DM particle mass-cw) plane, the mass of the DM 
particle is kept fixed and the value of the cross section 
for which 2.3 subhaloes are visible on average is searched. 
This corresponds to a 90% C.L. limit on av. 

Out of the subhaloes entering the histogram of Fig. [2] 
some are extended. Fig. [3] displays the distributions of 
fluxes and angular size of the subhaloes computed over 
100 random realizations of the Milky Way for two anni- 
hilation cross sections. Here the extension, dubbed 6*90, is 
defined as the angular size on the sky of the regions where 
90% of the annihilation signal arises from. Note that the 
brightest clumps are always more extended than the ave- 
rage, but they are never larger than about 1.5°. This is 
important because as it does not conflict with the on-off 
background subtraction technique discussed in Sec. IIIII 

To give an idea of how many of the subhaloes of Fig. [2] 
are visible for a given sensitivity, one can refer to the 
integrated luminosity function in Fig. 21 It presents the 
number of subhaloes in the field of view that have a flux 
larger than a given value versus the value of the flux. 
The error bars correspond to the rms of the distribu- 
tion due to the random nature of the halo realization. 
The results are presented for different particle DM confi- 
gurations : 500 GeV WIMP with annihilations into bb 
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Figure 3: Distributions of flux and angular sizes of subha- 
loes in the HESS field of view, obtained for annihilation cross 
sections of 10 -23 and 3 x 10 -26 cm 3 s -1 and stacked over 100 
random realizations of the Milky Way halo. 
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Figure 4: Integrated luminosity function for the VL-II sub- 
haloes in the HESS field of view. The considered WIMP mass 
is 500 GeV ; different cases for the annihilation parameters are 
considered : annihilations into bb and t + t~ and cross sections 
of 10" 23 and 10" 26 cm 3 s _1 . 



and t + t~ and cross sections of 10~ 23 and 10 -26 cm 3 s _1 . 
Should an experiment have an arbitrarily low flux sensi- 
tivity, it would lie on the lower-end asymptote and ob- 
serve all of the 168 objects on average. If the sensitivity 
map was constant over the whole field of view, to get 
a constraint one should just set the annihilation cross- 
section and mass of the DM particle to have 2.3 visible 
objects on average. In practice, the whole map including 
its fluctuations of sensitivity has to be used. 

Considering that HESS has an angular resolution of 
0.01°, a fraction (~ 50%) of the clumps that are found 
to be above the HESS sensitivity are extended. This has 



to be handled as the sensitivity map is built assuming 
point like sources. In that case, we make the conserva- 
tive choice to rescale the flux by lowering it to the value 
enclosed in the instrument angular resolution. Because 
the flux sensitivity concerns the integrated flux above 
some threshold, its value depends on an assumed spec- 
trum. Depending on the considered annihilation channel 
(hereafter \X ~ ► bb or t + t~) and on the DM particle 
mass, the values of the HESS sensitivity in each bin of 
the map are then properly rescaled. 

The left hand side of Fig. \E\ shows the 90% C.L. exclu- 
sion limit on <jv as a function of the DM particle mass. 
Two annihilation spectra arc considered : 100% bran- 
ching ratio annihilation channel in bb and t + t~, respec- 
tively, in order to somehow encompass all possible anni- 
hilation spectra for the DM particle. The limits on the 
annihilation cross section reach a few 10~ 24 cm 3 s _1 at 1 
TeV for the t + t _ spectrum. The dashed region corres- 
ponds to cosmologically relevant values for the annihila- 
tion cross section. The obtained constraints are 2 orders 
of magnitude above this region. 

These results are compared to those obtained with 
wide- field searches in Fermi data in [l2[ , which appear to 
be very complementary as they cover lower energies. For 
the sake of comparison, the best current constraints ob- 
tained with Cherenkov telescopes from targeted searches 
are from the Galactic center with HESS (of the order of 
10~ 24 - 10~ 23 cm 3 s _1 [l|) and Sagittarius dwarf galaxy 
(from 10 -25 to 10 -23 cm 3 s _1 , depending on the mode- 
ling of the source Q ) ■ The constraints obtained here from 
VL-II subhaloes are based on canonical assumptions and 
reach 10 -24 — IO -23 cm 3 s _1 ; they are thus among the 
most competitive to date obtained with 7 rays in this 
range of mass. 

In the case of an annihilation into gauge boson pairs 
(here T4 y+ VF _ ), it could happen that the cross section is 
significantly increased by means of the so-called Sommer- 
feld effect [H, [2^] . The enhancement factor S depends on 
the mass of the DM particle and on the relative velocity 
of the colliding particles, it ranges from a few percent up 
to very large values of ^10 4 . For V max <C c, the enhan- 
cement goes approximately as V~^ x before S reaches a 
plateau due the finite range of the Yukawa interaction. 
When m is close to a WIMP almost-bound state, the an- 
nihilation process becomes resonant with S oc V^ 2 ,,. The 
Sommerfeld effect has been modeled here in the case of 
an annihilation into W pairs, with an annihilation pro- 
ceeding towards the exchange of a 90 GeV boson and a 
coupling constant of g = 1/30. In the HESS sensitivity 
range, the resonances are obtained for m ~ 4.5 TeV and 
m ~ 17.6 TeV. In the specific case of annihilations within 
substructures, the enhancement can become very large, 
because the colder the subhalo, the larger the enhance- 
ment factor [3(| ■ In the subhaloes considered in VL-II, the 
maximal velocity Vmax ranges from 0.5 to 20 kms^ 1 . Ha- 
ving S oc V~^ x , the boost factor can indeed be large [3lj . 
In our Monte Carlo realizations, a specific Sommerfeld 
boost is assigned to each subhalo depending on its spe- 
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Figure 5: Left : Exclusion curves on av versus the DM particle mass m for HESS. The limit is calculated at the 90% C. L. for 
the DM clumps provided by the VL-II simulation. The DM particle is assumed to annihilate into purely bb and t + t~ pairs, 
respectively. The region of natural values of the velocity-weighted annihilation cross section of thermally produced WIMPs is 
also plotted. Right : Exclusion curves on av versus the DM particle mass m for HESS including the Sommerfeld enhancement 
effect. 



cific Vma X . The constraints obtained are displayed in the 
right panel of Fig.[3J Some predictions from supcrsymme- 
tric models with annihilation into W bosons, extracted 
from [32| [in the anomaly-mediated supersymmetry brea- 
king scenario (AMSB)], are also shown. These predictions 
do not include the S factor, so constraints on the unboos- 
ted cross section av/S are shown. Outside resonances, 
the limit is less than 2 orders of magnitude above the an- 
nihilation cross section expected for thermally produced 
WIMPS, but -thanks to the resonant Sommerfeld effect- 
a small region around 4.5 TeV is excluded. 



projection to the CTA are presented in Fig. [51 The ex- 
clusion limits are lower by a factor of ^10 than those 
obtained with HESS. In the conventional case (bb and 
t + t~, without Sommerfeld enhancement), they are rea- 
ching av values of a few 10 _25 cm 3 s~ 1 . In the case of 
Sommerfeld enhanced annihilations, some regions of the 
parameter space for the model could be excluded, since 
a large array of telescopes would have enough sensitivity 
to detect WIMPs in the mass range from ~3 to 6 TeV 
and close to the second resonance. 



V. PROSPECTS FOR CTA OBSERVATIONS 
PROGRAMS 

A. HESS-like Galactic plane survey 

The projected map for the CTA is used as in the pre- 
vious analysis of the HESS Galactic survey to make a 
prediction for the sensitivity of the future array. As a 
first step, the same field of view as HESS is used. We 
consider that a scan of the Galactic plane will for sure be 
performed by the CTA, so that this region of the sky 
is somehow the minimal guaranteed field of view. An 
exposure of 10 h in each pixels corresponds to a total 
observation time for building up the survey of 400 h. 
Concerning the extension of the sources, a slightly more 
optimistic method is used. Instead of rescaling the flux to 
the signal enclosed in the angular resolution, the whole 
signal is considered. Nevertheless, any clump that would 
be too much extended to allow for a proper background 
subtraction is excluded from the sample. Subhaloes with 
6>go > 1.5° are thus not considered. The results for the 



We conclude from Fig. [5J that using this field of view, 
the CTA will not be able to reach signals from the most 
natural WIMPs. One order of magnitude is gained with 
respect to HESS, but a factor of 2-10 is still necessary 
to reach the natural DM annihilation cross sections. An 
homogeneous increase of the exposure time will only im- 
prove the exclusion limits as the square root of exposure 
time in the background-limited regime, so one has to en- 
large the field of view instead of using longer exposure. In 
addition, the flux sensitivity along the Galactic plane will 
be limited by the population of newly detected sources at 
a flux level of 10 -12 cm _2 s~ 1 . The Galactic plane might 
also not be the best place to look for subhaloes since 
they could have been tidally affected by the disk. For 
those reasons, an observing strategy focusing on fields 
with absolute Galactic latitude of at least 0.5° should be 
preferred for DM subhalo searches, as it clearly appears 
in the lower panel of Fig. [TJ This is precisely the point 
developed in the next subsection. 
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Figure 6: Left : The projected sensitivity curves on av versus the DM particle mass m for a CTA Galactic plane survey. The 
limit is calculated at the 90% C. L. for the DM clumps provided by the VL-II simulation. The DM particle is assumed to 
annihilate into purely bb and t + t~ pairs respectively. The region of natural values of the velocity-weighted annihilation cross 
section of thermally produced WIMPS is also plotted. Right : Exclusion curves on av versus the DM particle mass m for a 
CTA-like array including the Sommerfeld enhancement effect. 
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Figure 7: Integrated luminosity function for the VL-II sub- 
haloes in the field of view of one quarter of the sky. The 
considered WIMP mass is 500 GeV ; different cases for the 
annihilation parameters are considered : annihilations into 6b 
and t + t~ and cross sections of 10 -23 and 10 -26 cm 3 s~ 1 



B. CTA survey of one-quarter of the sky 

To go further, one can note that larger scans of the sky 
will most likely be conducted by the CTA. In particular, 
a more extended survey of the order of a quarter-sky size 
is foreseen. In this section, the CTA sensitivity to DM 
annihilations is computed in the context of such an am- 
bitious program. A large survey increases the probability 
to find bright subhalocs in the field of view, which thus 
translates into better constraints. Such a survey should 
not include the Galactic plane where numerous sources 



are expected to shine and therefore decrease the sensiti- 
vity to DM clumps. On the other hand, the central region 
of the Milky Way is attractive since the VL-II subhalo 
distribution is peaked towards the center. For this study, 
the survey region is chosen to be from -90° to +90° in 
Galactic longitude and from -45° to +45° in Galactic lati- 
tude, excluding the Galactic plane between ±1.5°. Inside 
this region, the distribution of the number of subhalocs 
from the simulation is similar to the one of Fig. [2] but 
with an average value of 3907, and a rms of 324. The fact 
that this extended survey leads to better constraints is 
illustrated in Fig. [JJ where the same results as for Fig. @] 
are presented, but for the considered quarter-sky survey 
instead of the Galactic plane survey. Here again, the quo- 
ted error bars correspond to the RMS of the distribution 
over possible realizations. 

The sensitivity is taken to be constant on the en- 
tire field of view. Its value is calculated from the pre- 
vious CTA sensitivity map averaged for Galactic lati- 
tudes above 1.5° and corrected for a shorter exposure. 
A 5h exposure time in each pixel leads to a flux sensiti- 
vity of the order of 5x 10~ 13 cm -2 s -1 for a WIMP mass 
of 500 GeV. For this last study as well, the whole signal 
from extended subhalocs is considered, and objects lar- 
ger than 1.5° are rejected. As in previous sections, the 
value of the sensitivity is renormalizcd for each DM par- 
ticle mass. To a very good approximation, the decrease of 
the sensitivity due to the new population of extragalactic 
sources such as, e.g. active galactic nuclei, is negligible. 
The quoted value of the sensitivity is reached by poin- 
ting the whole array in the same direction, which is quite 
time consuming. Assuming a duty cycle of 1000 h of ob- 
servation per year, this quarter-of-the-sky survey can be 
completed within about 6 years of operation. An implicit 
assumption is that all unidentified sources do not present 
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Figure 8: Left : Exclusion curves on av versus the DM particle mass m for a CTA survey of a one fourth of the sky. The limit 
is calculated at 90% C. L. for the DM clumps provided by the VL-II simulation. The DM particle is assumed to annihilate 
into purely bb and t + t~ pairs respectively. The region of natural values of the velocity-weighted annihilation cross section 
of thermally produced WIMPs is also plotted. Right : Exclusion curves on av versus the DM particle mass m for the CTA 
including the Sommerfeld enhancement effect. 



the required features of DM source candidates. As for the 
Galactic survey, the observation of extension, variability, 
and energy spectrum will be used to rule out the uniden- 
tified sources as DM clump candidates. Notice however 
that should a couple of plausible DM sources be detected, 
additional dedicated observation time would be devoted 
for deeper studies. 

Figure E] shows the 90% C.L. exclusion limit on av as 
a function of the DM particle mass. Annihilation cross 
sections of a few 10~ 26 cm 3 s _1 are reached in the 200 
GeV - 3 TeV mass range in the case of annihilation into 
r pairs. In the scenario of an annihilation into W bosons 
pairs with Sommerfeld enhancement, all WIMPs from the 
AMSB with masses from 200 GeV to 6 TeV are within 
the reach of the CTA. 



VI. CONCLUSIONS 

We used for the first time a wide field-survey from Chc- 
renkov telescopes to constrain conventional DM substruc- 
ture scenarios. Unlike the case when DM annihilations 
are searched for towards selected sources, the constraints 
from this blind search do not rely crucially on the mode- 
ling of the DM distribution in the source. The constraints 



obtained out of the HESS Galactic plane survey are still 
2 orders of magnitude higher than the thermal WIMP 
region. Thus most natural models for WIMPs as dark 
matter are out of the reach of current generation ground- 
based Cherenkov telescopes with wide-field surveys and 
realistic observing time. However, the limits reached in 
the av — m plane are very competitive compared to other 
strategies such as targeted searches. By using the same 
Galactic plane field of view, we show that the discovery 
of particle DM in the form of WIMPs is unlikely to be ac- 
cessible for the next generation of Cherenkov telescopes 
such as the CTA. However, by considering an ambitious 
but realistic quarter-of-the-sky survey with the CTA, it 
is shown that the thermal WIMP promised land can be 
hit. Note that such a survey will likely be conducted by 
CTA independently of particle DM considerations. This 
search for DM subhaloes will therefore not be in conflict 
with other physics programs. 
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